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Abstract. The only invariant speed in special relativity is c; therefore, if some neutrinos 
travel at even tiny speeds above c, normal special relativity is incomplete and any super- 
luminal speed may be possible. I derive a limit on superluminal neutrino speeds v ^ c at 
high energies by noting that such speeds would increase the size of the neutrino horizon. 
The increased volume of the Universe visible leads to a brighter astrophysical neutrino back- 
ground. The nondetection of "guaranteed" neutrino backgrounds from star- forming galaxies 
and ultrahigh energy cosmic rays (UHECRs) constrains v/c at TeV-ZeV energies. I find that 
u/c < 820 at 60 TeV from the nondetection of neutrinos from star- forming galaxies. The 
nondetection of neutrinos from UHECRs constrains v/c to be less than 2500 at 0.1 EeV in a 
pessimistic model and less than 4.6 at 4 EeV in an optimistic model. The UHECR neutrino 
background nondetection is strongly inconsistent with a naive quadratic extrapolation of the 
OPERA results to EeV energies. The limits apply subject to some caveats, particularly 
that the expected pionic neutrino backgrounds exist and that neutrinos travel faster than c 
when they pass the detector. They could be improved substantially as the expected neutrino 
backgrounds are better understood and with new experimental neutrino background limits. 
I also point out that extremely subluminal speeds would result in a much smaller neutrino 
background intensity than expected. 
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1 Introduction 

An early question with implications for cosmology is Olber's Paradox: Why is the sky dark 
at night? In a static, infinite Universe every line of sight should eventually intercept the 
surface of a star, and so the sky should be as bright as a stellar surface (see [1] for a critical 
review of the early history of Olber's paradox). The resolution of the paradox arises from the 
finite age of the Universe in Big Bang cosmology: stars only radiate for finite times, limiting 
the energy density of the background; in terms of flux, the Universe has a finite age and so 
we see only those galaxies where the light has had time to reach us [2].^ The background flux 
from the sky therefore depends on the speed v of the particle, which sets the length scale of 
the "horizon" for that particle. 

The recent claim of superluminal neutrinos {Av/c = t;/c — l~2x 10~^ at ~ 17 GeV) 
from the OPERA experiment [4] has excited much interest. Superluminal particles have 
been predicted to arise through Lorentz invariance violation or propagation through extra 
dimensions. Several constraints on superluminal particles exist. The most powerful limits 
use the kinematics of superluminal particles to infer new energy loss processes which would 
limit the range of particles, or alter known particle reactions. Such limits have previously 
been derived for protons [5] and photons [6] . Recent limits of this kind for neutrinos include 
those from e^e~ pair production [7~10], neutrino splitting [11], and suppression of pion decay 
[10, 12]. These limits rule out the OPERA results in standard Lorentz invariance theories, 
and they become stronger at high energy. But it is conceivable that these limits can be 
evaded, for example, if photons and electrons have the same speed v ^ c at high energies as 
neutrinos [8] , if neutrinos transform into a superluminal particle that does not couple with 
the Standard Model sector [13] , or in some proposed modifications of Lorentz invariance that 
alter the laws of energy- momentum conservation [14] or the spacetime metric [15]. 

Generic limits arise from timing considerations if particles are detected from an astro- 
physical transient [16]. In particular, gamma-ray bursts last only a few seconds in our frame 
and are billions of light years away; even a small difference in propagation speed will accumu- 
late into huge differences in arrival times. Fermi-LAT observations of gamma-ray bursts has 
already severely constrained velocities of the form Av/c oc E for photons from the detection 
of GeV 7-rays coincident with 7-ray bursts [17]. 

^In addition, while the Universe once was as bright as a stellar surface, the cosmic expansion redshifted 
that radiation to microwave frequencies (e.g., [3]). However, the cosmic expansion plays a relatively minor 
role in setting the energy density or energy flux of starlight [2]. 
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At TeV energies, the Universe becomes opaque to 7-rays, but neutrinos can also provide 
powerful constraints if they are ever detected [18, 19]. Unfortunately, no astrophysical neu- 
trinos have been detected except from the supernova SN1987A and the Sun. The detection of 
core collapse neutrinos from SN1987A a few hours before it became visible in light provides 
strong constraints on vjc at tens of MeV [20]. 

In special relativity, any speed minimally faster than light in one frame can appear 
much faster than light in another frame, suggesting that once we allow the possibility of 
superluminal particles at all, we should consider particles much faster than c. Indeed, some 
ideas for superluminal neutrinos predict that the neutrino speed should grow at higher energy. 
Higher order effects may cause A?;/c to asymptote to some small limit, but once we allow the 
possibility of superluminal particles, there may be nothing special about the photon speed c 
per se so that high energy neutrinos travel at speeds much faster than c. For all we know, 
higher order terms may cause Af/c to grow even faster at some special energy scale (e.g., 
exponentially) and asymptote to u 3> c. This suggests we should look at higher energies 
where superluminal speeds could be huge. 

Furthermore, if neutrino speed limits from SN1987A and the results of the OPERA 
experiment apply to all flavors, then Ai)/c must grow quickly at least between the MeV 
scales of SN 1987A's neutrinos and the GeV scales of the OPERA neutrinos [21]. Assuming 
Au/c oc the OPERA result implies Au/c !^ 1 at 3 TeV; if l^vjc oc E, then the OPERA 
result implies ^vjc ~ 1 at 0.7 PeV. But regardless of whether the OPERA experiment result 
is correct, any ISxijc that can be directly measured at low energies may be extrapolated to 
u ^ c at high energies with a strong enough energy dependence. 

I present a constraint on extreme superluminal neutrino speeds at very high energies, 
under certain assumptions described in the next section. In short, the faster neutrinos can 
travel, the more of the Universe we can see with neutrinos. Thus the flux intensity of the 
diffuse background of all neutrino sources in the Universe should increase linearly with the 
neutrino speed in accordance with Olber's Paradox. However, there already are limits on 
the diffuse neutrino background, and some of these are already near theoretical predictions for 
the neutrino background assuming v = c. The fact that we do not see a neutrino background 
at high energies therefore implies that v is not many orders of magnitude higher than c at 
these energies. 

2 Explanation and Assumptions 

The idea behind this limit is that superluminal neutrinos increase the intensity (defined here 
as number flux per steradian) of the diffuse background. Suppose neutrino sources have equal 
luminosity and uniform number density throughout the Universe. In a shell of thickness dR 
at distance R from Earth, the flux of neutrinos per source goes as R~^ but the number of 
sources in the shell goes as R^dR. Thus each shell has the same flux as viewed from Earth, 
and the integrated flux increases linearly with the maximum observable distance. Neutrinos 
with V ^ c would be observable out to a larger horizon, and the neutrino background would 
have a correspondingly higher intensity. 

Note the number density of the neutrino background is not increased: that is set simply 
by the rate of neutrino injection per comoving volume (ultimately from the energetics of 
neutrino production) and the age of the Universe. Rather, the higher speed of neutrinos 
increases the rate at which the neutrinos hit a detector on the Earth. Thus, models in which 
the neutrinos transform into some intermediate state which is superluminal and then turn 
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back into subluminal particles before hitting the detector (e.g., [13]) are not constrained — 
those models alter which sources contribute to the background but not its intensity. 

The number density spectrum dn/dE [cm~^GeV~^] of neutrinos is calculated as / dtpcom[dQ / dE'][dE' / c 
where dt = dz/[HQ{l + z)-\/Q.j\^ + ^m{1 + z^] describes the lookback time, pcom is the co- 
moving density of sources, and dQ/dE' is the source-frame number luminosity [s~^GeV^^] 
of a source. I have assumed that neutrinos accumulate over cosmic time in the FRW metric, 
rather than some other coordinate time. From the number density spectrum, we can calcu- 
late the present-day background intensity simply as d^/dE = dn/dE x v/{Att), where v is 
the current neutrino speed at energy E, to get 



I will assume that dE' /dE = (1 -|- z)\ this holds for particles with E k, pc{\ + 5), 
which leads to dE/dp = v ~ c{l + 6). The redshift effects are relatively minor unless the 
energy dependence is much more extreme than E' (x {1 + z) [2]. It can then be shown that 
d^/dE scales directly with the constant superluminal velocity at all energies: thus as v — > oo, 
d^/dE — )■ oo in accordance with Olber's Paradox. 

However, some caveats apply to these limits: 

1. As discussed above, I assume that neutrinos themselves are superluminal in the space 
the detector occupies, because this is a limit on flux and not density. 

2. I assume that expected astrophysical neutrino backgrounds from pion production pro- 
cesses actually exist; if pions do not decay or neutrinos decay en route [7, 10, 12, 19], 
the limits will not apply. 

3. I assume that the neutrino speed is constant at every point in the Universe with respect 
to the frame in which the cosmic expansion is isotropic. If this is not the case, there 
is a preferred location in the Universe. Elsewhere, the diffuse neutrino background 
will appear anisotropic and superluminal neutrinos may come from our future with 
respect to cosmic time (from — 1 < z < 0). Furthermore, the neutrino density would 
represent the integral of energy injection over some coordinate time that may bear 
little resemblance to cosmic time, and the neutrino density could then be far lower 
than expected. 

4. Likewise, I ignore the effects of the Earth's motion to this frame, which will introduce 
large anisotropy in the Earth-frame speed iiv/c> c/v,^ ^ 800 for = 370 km s~^ [22]. 
The Earth-frame number density will also appear to be anisotropic, with dn/dO — t- 
as v{6)/c —7- oo (because neutrinos arrive the moment they are emitted in Earth-frame 
and there are only a finite number of neutrinos emitted within the neutrino horizon) 
so the measured event rate is not changed to first order. 

5. I assume that the region of the Universe with neutrino sources is infinitely large: if it 
is not, the enhancement of the background will reflect the true size of the neutrino- 
emitting Universe and not the neutrino horizon (although this may be interesting in 
its own right). 




pcom 



dQ dE' 
dE'liE 



(2.1) 
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6. Finally, I assume the Universe is infinitely transparent to neutrinos. However, at en- 
ergies > 10^^ eV, neutrinos may be absorbed over distances much greater than c/Hq 
through interactions with low energy neutrinos and photons [23].^ 

3 Application 

While we have not detected TeV or higher energy neutrinos from astrophysical sources, we 
expect that neutrino sources exist because of the observations of cosmic ray (CR) protons and 
nuclei. These particles create pions through either collisions with ambient nucleons {pp) or 
photons (p7). Then, assuming the kinematics are not altered, the pions decay into secondary 
electrons and positrons, 7-rays, and neutrinos. Thus wherever there are CRs and enough 
gas or radiation (of sufficient energy for the p^ process to occur) , we expect neutrinos to be 
emitted. 

I consider GeV-PeV neutrinos from star-forming galaxies and PeV-ZeV neutrinos from 
extragalactic UHECRs, which both probably exist, though the exact levels of these back- 
grounds is still unclear. In addition, there may be intense neutrino backgrounds from active 
galactic nuclei [24] and other sources, but the existence of these backgrounds is more uncer- 
tain. Finally, I note that the MeV diffuse supernova neutrino background is expected to have 
an intensity only a factor of ~ 4 below current limits. Thus f/c < 4 at ~ 10 — 30 MeV [25], 
although these limits are unfavorable compared to extant ones from SN1987A [20] and the 
OPERA experiment at GeV energies [4]. 

3.1 The Star-Forming Galaxy Neutrino Background 

The Milky Way glows in pionic 7-ray emission, indicating the presence of CR protons in- 
teracting with gas throughout the Galactic disk [26]. The detection of synchrotron radio 
emission [27] and now GeV 7-ray emission in other nearby galaxies [28, 29] indicates that 
CRs are also accelerated in them as well, and furthermore, that the CR energy injection rate 
is related to the star-formation rate. This led to the prediction of a minimum "guaranteed 
7-ray background" from pp collisions throughout galaxies in the Universe [30, 31], but the 
same process also produces a neutrino background [32, 33]. At high energies the neutrino 
spectrum of a galaxy is proportional to, and of roughly the same intensity as its 7-ray spec- 
trum [32]. I consider two subclasses of star- forming galaxies: "normal" (Milky Way-like) 
galaxies and starburst galaxies. 

In normal galaxies, CR protons escape by diffusion, leading to rapidly falling steady- 
state CR proton, 7-ray, and neutrino spectra {dN/dE oc E^'^''^; [34]). The spectrum steepens 
at a few PeV (the "knee"); since the typical energy of a neutrino is about 5% of the proton 
energy, this translates to a cutoff in the neutrino spectrum at ~ 100 TeV [33]. I therefore 
assume a E~'^''^ power law spectrum for the Milky Way, with a cutoff at 100 TeV. Ref. [35] 
modeled CR propagation in the Milky Way and found a > 1 GeV pionic 7-ray luminosity of 
3.7 X 10^^ ergs s~^. Scaling to the star- formation rate of the Milky Way SFRmw = 2 Mq yr"^ 
(e.g., [36]), I use a neutrino spectrum of 

g.....03V.C.V- (4)-"e.,(^)(^). 

■^A similar absorption-limited background was proposed by Olber as the explanation for why the night sky 
is dark in optical light. 
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Table 1. Summary of constraints on vjc of neutrinos of energies from 30 GeV to 300 TcV from 
star-formation neutrino background. At 100 TeV and above, the neutrino background depends on 
where the "knee" in the CR nuclei spectrum is. I conservatively compare the spectra to the atmo- 
spheric neutrino background, and also compare the starburst component alone to the E~^'^ limit from 
AMANDA. 

Limit on v/c (Atmospheric) Limit on v/c AMANDA) 



Energy 


Normal Galaxies 


Starbursts 


Total 


Starbursts 


30 GeV 


7.6 X 10^ 


5.3 X 10^ 


3.1 X 10^ 




100 GeV 


6.5 X 10^ 


2.5 X 10^ 


1.8 X 10^ 




300 GeV 


4.3 X 10^ 


9.5 X lO'' 


7.8 X 10^ 




1 TeV 


2.7 X 10^ 


3.2 X 10"^ 


2.8 X 10^ 




3 TeV 


1.5 X 10^ 


1.1 X 10^ 


1.0 X 10^ 




10 TeV 


7.4 X 10^ 


2900 


2800 


100 


30 TeV 


4.4 X 10^ 


990 


970 


150 


100 TeV 


7.4 X 10^ 


980 


960 


550 


300 TeV 


6.5 X 10^ 


5500 


5400 


1.1 X lO'' 



Starburst galaxies have recently been detected in GeV-TeV gamma-rays [29] , and these 
gamma-rays are believed to be pionic [37, 38]. The high star- formation rate in the starburst 
regions leads to high cosmic ray energy densities; when combined with the high gas densities 
in the starbursts, high pionic luminosities are expected [31, 39]. I assume a spectrum for 

starburst galaxies as observed, and scale to the > 1 GeV luminosity of M82 (1 .9 x lO^'' ergs s"^ ; 
[38]) and M82's star-formation rate SFRm82 = 8 Mq yr^"*^ from its infrared luminosity [40] ^: 

^ = 7.6 X lO-s-'GeV- f^V'.xp ( ( J^] . (3.2) 

dE' \GeV J ^\100 TeV J \Mq yr-i J ^ ' 

The cosmic star-formation rate is given in ref. [41]. In order to calculate the spectrum of 
the neutrino background, the fraction of star-formation in normal galaxies versus starbursts 
is necessary. While debated in the literature [42], I use a conservative 5% at all redshifts in 
starburst galaxies [43], with the rest in the fainter normal galaxies. Note that I ignore the 
effects of gas evolution, which can enhance the neutrino flux in normal galaxies [33, 44]. I 
assume Q.m = 0.3, = 0.7, and Hq = 70 km s~^ Mpc~^ for this calculation. 

Current neutrino limits with IceCube are for an E~'^ spectrum and only apply at en- 
ergies of 35 TeV — 7 PeV [45], where the predicted neutrino spectrum tails off (Figure 1). 
However, the star- formation neutrino background is certainly smaller than the atmospheric 
neutrino background, which is measured up to 400 TeV [46]. I therefore take the ratio of the 
atmospheric neutrino spectrum from [47] and the predicted neutrino background to find a 
constraint on v/c. This is very conservative, since pionic astrophysical neutrino backgrounds 
have different spectra, flavor ratios, and angular distributions than atmospheric neutrinos. 

These constraints on v/c are weak, especially at low energies where the atmospheric 
background is high (Figure 1 and Table 1). At 100 GeV, v/c < 2 x 10^, but the best 
constraint is at 63 TeV, where v/c < 820. The constraints would be stronger if the starburst 
fraction were higher (the starburst contribution already dominates above 12 GeV); if it were 
100%, the best constraint would be f/c < 42 at 63 TeV (compared to ~ 340 from a naive 
quadratic extrapolation of the OPERA result). 

^As converted into a "Salpeter A" initial mass function of stars, as used in [41] 
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[TeV] 

Figure 1. Predicted pionic neutrino fluxes from star- forming galaxies with v/c = 1 (black) compared 
to the atmospheric neutrino spectrum (thick black; [47]) and the AMANDA limits on a E~^'^ spectrum 
from 6.3 TeV to 1.3 PeV (thick black; [48]). The dotted line is the normal galaxy contribution and 
the dashed line is the starburst contribution (if 5% of the cosmic star formation rate is in starbursts) . 
The grey line is the background scaled up by the maximum allowed v/c — 820. 

However, as IceCube sets increasingly stringent limits on neutrino backgrounds at 

10 TeV, the limits on v/c should decrease as the starburst neutrino background approaches 
detection. AMANDA set a limit on E"'^-'^ spectra extending from 6.3 TeV to 1.3 PeV [48]. 
Comparing to my predicted starburst spectrum and ignoring the fact that its energy cutoff 
would alter the neutrino analysis, I find v/c < 90 at 6.3 TeV. 

3.2 The UHECR Neutrino Background 

Ultrahigh energy cosmic rays (UHECR) are observed at energies up to lO^'' eV and beyond, 
and are generally believed to be extragalactic in origin, although their composition is un- 
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Table 2. Summary of constraints on vjc of neutrinos of energies from 100 GeV to 100 TeV from 
UHECR neutrino background. The GZK models are from ref. [58]: the optimistic model has a proton 
composition and strong evolution; the medium case has a proton composition and weaker evolution; 
the pessimistic case has an iron composition. Also included are the bounds if the neutrino injection 
rate is at the Waxman-Bahcall limit [59]. 

Limit on i)/c 







GZK Models 




Waxman-Bahcall 




Optimistic 


Medium 


Pessimistic 


Non-Evolving 


With Evo 


3 PeV 


75 


170 


6.2 X 10^ 


15 


5.0 




1 on 


L 1 U 




23 


7.8 


30 PeV 


220 


570 


6.9 X lO'' 


35 


12 


100 PeV 


130 


570 


2.7 X lO'' 


27 


9.1 


300 PeV 


28 


190 


9700 


21 


7.1 


1 EeV 


6.6 


52 


3100 


12 


3.9 


3 EeV 


4.9 


35 


2800 


9.4 


3.1 


10 EeV 


5.2 


56 


1.1 X lO'' 


8.6 


2.9 


30 EeV 


8.7 


220 


2.9 X 10^ 


10 


3.4 


100 EeV 


13 


1500 




7.0 


2.3 


300 EeV 


82 


6.5 X 10^ 




8.2 


2.7 


1 ZeV 








11 


3.6 


10 ZeV 








25 


8.3 


100 ZeV 








88 


29 



known. As they traverse the intergalactic medium, UHECR nucleons interact with the CMB 
to produce pions above the GZK (Greisen-Zatsepin-Kuzmin) threshold eV; [49]); 

these presumably ultimately decay into neutrinos and other particles [50]. The existence of 
the photopion process is supported by the suppression of the UHECR flux above 40 EeV 
[51], although photodisintegration of heavy ions could also create this dip [52]. 

A number of experiments have sought these GZK neutrinos, and stringent limits have 
been placed, ruling out some models already (e.g., [53-57]). Even in the most pessimistic 
models where UHECRs consist of entirely heavy nuclei, some neutrinos are expected from 
photopion production on the extragalactic background light in the infrared, as well as pho- 
topion production from secondary protons after photodisintegration of heavy nuclei [58]. A 
large vjc at EeV energies would enhance the flux from these small backgrounds to the point 
where they are detectable. 

I consider three models of the UHE neutrino flux from [58] , who consider a wide range 
of parameters. In their most optimistic case, UHECRs are protons and injection evolves 
strongly with redshift (as radio galaxies do), in a more typical case, UHECRs are protons 
and evolve more weakly with redshift (as star formation does), and in a pessimistic case, 
UHECRs are iron nuclei and injection does not evolve with redshift (Figure 2). 

Using these predictions, I find that the best constraint of all energies is ?;/c < 2500 at 
0.18 EeV in the pessimistic model, t;/c < 35 at 3.2 EeV in the typical model, and v/c ^ 4.6 
at 4.2 EeV in the optimistic model. These constraints are particularly interesting in that a 
naive quadratic extrapolation of the OPERA result would indicate v/c ~ 9x 10^^(i?/l EeV)^; 
even a linear extrapolation gives u/c ~ 1500(-E/1 EeV). 

Finally, a characteristic density for the UHE neutrino background is the Waxman- 
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[10'9 eV] 

Figure 2. Predicted neutrino fluxes from UHECRs with v/c = 1 (black) compared to instrument 
sensitivities (thick solid black). The dotted line is a pessimistic pure iron composition model, the 
short-dashed line is a typical proton dip model, and the long-dashed line is a proton model with 
strong evolution from Ref. [58]. The dash-dotted lines are the Waxman-Bahcall bounds on the 
neutrino background [59]. In grey, I show the models scaled up by the maximum v/c allowed by the 
data. Observed limits are from AMANDA [53, 54], IceCubc [54], Auger [55], RICE [56], and ANITA 
[57]. In grey, I also show the projected limits from ARA [62] and SKA [63]. 

Bahcall bound, in which the energy injection rate in UHE neutrinos is less than that of 
observed UHECRs [59]. If the UHE neutrino energy injection rate is at the Waxman-Bahcall 
limit, then v/c\s less than 90 (if there is no redshift evolution) or 30 (with redshift evolution) 
over the entire range from 2 PeV to 100 ZeV. The most stringent limit would then be at 
140 EeV: v/c < 6.8 without evolution or v/c < 2.3 with evolution. However, the neutrino 
background density could easily be much lower than the Waxman-Bahcall limit, which would 
weaken these limits. 
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4 Conclusion 



I have argued that if neutrino are superluminal and stable, it is conceivable that f ^ c at 
TeV energies and beyond. In standard Lorentz invariance violation models, this possibility 
is ruled out [7-11], but some have suggested that there could be ways around these limits 
[13-15]. Astrophysical neutrinos could provide strong tests of superluminal neutrino speeds, 
especially through the timing of neutrinos arriving from transients, but no astrophysical 
neutrinos of these energies have been detected yet. 

By taking advantage of the increased neutrino horizon size if ^ c, I constrain v/c on 
TeV-ZeV scales. The lack of neutrinos from star- forming galaxies conservatively constrains 
v/c < 820 at 60 TeV and the lack of GZK neutrinos conservatively limits f/c < 2500 at 0.1 
EeV energies. The limits are subject to several caveats: most notably, pionic neutrinos must 
actually be produced and survive, and must actually be superluminal when they cross the 
detector since the limit is on background intensity and not density. 

The limits can be improved by further understanding of neutrino sources, so that we 
have a more accurate baseline prediction to compare the limits to. A determination of the 
composition of UHECRs [60] and constraints from the cascade radiation of GZK photons 
(e.g., [61]) will allow us to better predict what the "guaranteed" GZK neutrino flux is. 
Likewise better understanding of the neutrino background from star-forming galaxies may 
come from studies of the accompanying pionic gamma-ray background. 

On the observational side, experiments such as IceCube, the Askaryan Radio Array 
(ARA; [62]), and the Square Kilometer Array (SKA) will push down the limits on the neutrino 
background, with ARA expected to detect the UHE neutrino background except in the 
most pessimistic cases. Future neutrino detections of sources like AGNs may demonstrate 
additional "guaranteed" neutrino sources. Eventually, if the uncertainties in the predicted 
neutrino background are less than a factor of 2 and the neutrino background is detected at 
predicted levels, we will be able to probe directly the regime where Av ~ c where any change 
in energy dependence might naturally occur. 

Of course, the general idea behind this method is not just limited to neutrinos, but 
could apply to any particle. On the other hand, just as high order effects might cause v 
to increase rapidly at high energies, by the same logic, high energy neutrinos may travel 
far slower than c. Extreme subluminal speeds at high energies would have opposite effects 
as superluminal speeds: the background flux (but not density) would be much lower than 
predicted. A smaller horizon is already predicted for the relic cosmic neutrino background, 
where the neutrino mass is important [64]. 

If it turns out that v really is much larger than c, a number of interesting consequences 
arise in neutrino astronomy. We would see photon-observed galaxies as they are now in cosmic 
time, millions or even billions of years after their light was emitted, in completely different 
stages of their evolution. Since there is less star formation [41] and active galactic nucleus 
activity [65] at z ~ 0, we would see far fewer active galaxy neutrino sources and gamma 
ray burst neutrino flares within the photon horizon, even in a blind survey not affected by 
long delays. The diffuse neutrino background is brighter only because of the contributions 
of neutrinos emitted at z > 1 from sources far beyond the photon horizon that individually 
have much lower observed fluxes. In contrast, if f ^ c, then we would see sources near the 
neutrino horizon as they were at high redshift, when they presumably were brighter, even 
though the background whole is smaller. 
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